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ABSTRACT 





The method of distributed calibration of a probe microscope scanner consists in a search for a net of local 
calibration coefficients (LCCs) in the process of automatic measurement of a standard surface, whereby 
each point of the movement space of the scanner can be defined by a unique set of scale factors. Feature- 
oriented scanning (FOS) methodology is used to implement the distributed calibration, which permits 
to exclude in situ the negative influence of thermal drift, creep and hysteresis on the obtained results. 
The sensitivity of LCCs to errors in determination of position coordinates of surface features forming 
the local calibration structure (LCS) is eliminated by performing multiple repeated measurements fol- 
lowed by building regression surfaces. There are no principle restrictions on the number of repeated LCS 
measurements. Possessing the calibration database enables correcting in one procedure all the spatial 
distortions caused by nonlinearity, nonorthogonality and spurious crosstalk couplings of the microscope 
scanner piezomanipulators. To provide high precision of spatial measurements in nanometer range, the 
calibration is carried out using natural standards — constants of crystal lattice. The method may be used 
with any scanning probe instrument. 
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1. Introduction 


Usually, a probe microscope scanner is characterized by three 
calibration coefficients Ky, Ky, Kz representing sensitivities of X, 
Y, Z piezomanipulators, respectively (to take into consideration 
a possible nonorthogonality of X, Y piezomanipulators, an obliq- 
uity angle should be additionally determined) | 1-4]. Because of 
piezomanipulator’s nonlinearity [5-7] and spurious crosstalk cou- 
plings, the probe microscope scanner may be described by the 
above coefficients only near the origin of coordinates, where the 
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influence of the distortion factors is insignificant. As moving away 
from the origin of coordinates, the topography measurement error 
would noticeably increase reaching the utmost value at the edge of 
the scanner field [8]. 

The problem may be solved by using a distributed calibra- 
tion, which implies determining three local calibration coefficients 
(LCCs) Kx, Ky, Kz for each point of the scanner movement space, 
which can be thought of as scale factors for axes x, y and Z, respec- 
tively [9-13]. Areference surface used for calibration should consist 
of elements, called hereinafter features, such that the distances 
between them or their sizes are known with a high precision. 
The corrected coordinate of a point on the distorted image of an 
unknown surface is obtained by summing up the LCCs related to 
the points of the movement trajectory of the scanner (see Sections 
2.2 and 3). 
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Both lumped and distributed calibration of the probe micro- 
scope scanner should be carried out by the data where distortions 
caused by drifts (thermal drifts of instrument components plus 
creeps of piezomanipulators) are eliminated. Otherwise, the mea- 
surements will have large errors | 1—3,14—22]. In the present work, 
to eliminate the negative influence of thermal drift and creep on 
the distributed calibration results, the methods are used of feature- 
oriented scanning (FOS) and of counter-scanning suggested in Refs. 
| 17-20] (see Section 2). 

It is necessary to note that manual calibration is quite a labori- 
ous process even in case the scanner is characterized by the three 
parameters. As to the search for distributed coefficients, it is only 
possible under full automatic control. 

The investigation of the distributed calibration technique as 
a whole consists of three parts. In the first part presented here, 
the task of distributed calibration is formulated and its solution 
based on the FOS approach is described. Simultaneously with the 
description of the performed operations, key notions and tech- 
niques underlying FOS are shortly introduced and explained. The 
second part of the investigation is given in Ref. [14], it is devoted 
to the so-called virtual calibration mode. With this mode, instead 
of measurement of the real surface of a standard, the calibration 
program performs “measurement” of an image of the standard sur- 
face, which has been obtained earlier during a regular scanning. The 
virtual mode is intended for simulating the process of calibration 
and validating the analytical solutions found in the present work. 
The use of the virtual mode allows for analyzing the operation of 
a probe microscope scanner, detecting and estimating the errors 
that occurred. The third part of the investigation is given in Ref. [8], 
where the experimental results obtained during real distributed 
scanner calibration by crystal lattice of highly oriented pyrolytic 
graphite (HOPG) are presented and discussed. 


2. Measurement process of a standard surface 


Shortly, the suggested method of distributed calibration [11,12] 
of a probe microscope scanner is as follows. First, the scanning field 
is “covered” with a square cell net where the nodes correspond to 
the absolute integer coordinates x, y, z of the scanner (see Fig. 1). 
During the calibration, a determination of LCCs is carried out in the 
neighborhood of each net node. For this purpose, a sequence of the 
three basic measurement operations is carried out at a standard 
surface: probe attachment, aperture scanning, and skipping (see 
Fig. 2, pos. 2-4) [17]. 

The attachment of the microscope probe to a surface feature 
allows capturing the feature located near the current net node and 
to holding it within the “field of view” of the instrument. Scanning a 
small neighborhood called an aperture around the captured feature 
followed by recognition of topography in the obtained scan enables 
detecting local calibration structure (LCS). LCS usually consists of 
three features A, B and C such that the distances between them are 
a priori known precisely. To get the sought-for LCCs, it is required 
to measure the distances between the features A and B, Aand C [1]. 
Drift-insensitive measurement of these distances is carried out by 
means of the skipping operation. 


2.1. Detailed description of the calibration operations 


Let us consider the process of the distributed calibration in 
detail. The microscope probe moves by the net nodes like by a raster 
(see Fig. 1(a)). The size of the net square side approximately defines 
how smoothly the LCCs are changing within the found distribution. 
A pause is inserted after each movement (see Fig. 2, pos. 1), dur- 
ing which the piezoscanner creep produced by this movement is 
decaying. The larger the distance between the nodes of the initial 
net is, the longer the pause is set. 











Fig. 1. Partitioning of the scanner movement space with a net, LCCs being sought 
for in the vicinity of each node. The arrows show a spatial movement trajectory by 
the net nodes during the distributed calibration. (a) Calibration in the lateral planes 
with “fixed” positions of the scanner Z manipulator. (b) Calibration in the vertical 
planes with fixed positions of the scanner Y manipulator. In both cases, position of 
the scanner Z manipulator is set up by means of a coarse approach stage. First, a local 
scanning of the standard surface is executed in the vicinity of each net node, then 
the nearest LCS is detected, by which the LCCs are determined. In order to reduce 
the resulting creep, the movement trajectory from one node to the other is chosen 
so that the movement in the adjacent (a) lines, (b) columns and in the adjacent (a) 
horizontal, (b) vertical planes be implemented in the mutually opposite directions. 


2.1.1. Probe attachment 

When the pause is over, the feature A of the standard sur- 
face nearest to the current net node (x, y) is being “captured” (see 
Fig. 3(a)). Topography elements that look like hills or pits may be 
used as features of the standard surface. Those may be for example: 
atoms, interstices, molecules, clusters, grains, nanoparticles, crys- 
tallites, quantum dots, nanoislets, nanopillars, pores, etc. [17-21]. 

The feature capture (tracking) is carried out by scanning a small 
square neighborhood around the feature, recognizing surface fea- 
tures | 1,17] on the obtained local scan and then moving the probe 
to the position of the feature located nearest of all to the local scan 
center. The described sequence of operations is known as probe 
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3 Aperture scanning and recognition 
of local calibration structure (LCS) 


4 Skippings between LCS features 


SS) Calculation of local 
calibration coefficients (LCCs) 


6 Store results in 
calibration database (CDB) 


t Move X/Z manipulator to 
next absolute position x/z 
8 Move Y/X manipulator to 
next absolute position y/x 
g Move Z/Y manipulator to 
next absolute position Z/y 
10 Repeat / times 


Fig. 2. Simplified flowchart of distributed calibration of the probe microscope scan- 
ner. The designations in pos. 7-9 at the left of the slash correspond to calibration 
in the lateral planes with “fixed” positions of Z manipulator (see Fig. 1(a)); the ones 
at the right correspond to calibration in the vertical planes with fixed positions 
of Y manipulator (see Fig. 1(b)). Absolute position of the scanner Z manipulator is 
adjusted by means of a coarse approach stage. LCCs are searched for in the vicinity of 
absolute integer coordinates x, y, z of the scanner. In some practical cases, pause and 
probe attachment may be omitted. To increase accuracy, the calibration is repeated 
! times. The number | of repeated calibrations has no limitations. In the CDB stored 
are: the absolute real coordinates of points of the scanner space for which the LCCs 
were obtained and the LCCs values themselves corresponding to these points. 







attachment (see Fig. 2, pos. 2) [17]. The captured feature is then 
being held for some time within the field of view of the instrument 
by means of successively repeated attachments. 

When the feature is being captured for the first time, a square 
area is scanned out of such size that being located arbitrary relative 
to the structure of the standard surface it is able to enclose at least 
one feature of the standard surface (see Fig. 3(a), “1st Attachment”). 
To increase the productivity, while executing next attachments, the 
sizes of the scan area is reduced to the sizes of a segment, which is 
a square scan enclosing one feature only (see Fig. 3(a)) [17]. 

During the attachments, the drift velocity is determined 
[17-21]. If the drift velocity turns out to be greater by module than 
a certain preset value then the attachment is repeated. Usually, 
the preset value is chosen to represent the mean velocity typical 
of the given microscope with which the microscope drifts after 
warming up [8,14]. Thus, the attachment allows to determine the 
end moment of the creep induced by movement to the next net 
node. In case of a large net step and corresponding strong creep 
excitation, it is recommended to set a longer pause in order to 
prevent substantial deviation of the scanner from the accepted 
movement trajectory by net nodes (see Fig. 1) while performing 
series of attachments. The fact is that with a substantial deviation, 
the following scanner return to this trajectory will excite a strong 
creep in its turn. 


(a) 


. 


ment: 





Fig. 3. Schematic image of the standard surface. Nodes of the net dividing the scan- 
ner field into square cells are designated as “+”. The nodes correspond to the absolute 
integer coordinates of the scanner. The point (x, y) is the current net node, the LCCs 
are searched within its vicinity. The group of three features A, B and C of the standard 
surface serves as the LCS, the distances between the features or sizes of the features 
are known a priori with a high accuracy. Positions of the features detected during 
calibration are designated as “*”. LCCs are determined for the “gravity center” of the 
LCS ABC designated as “e”. The attachment procedure scans a square neighborhood 
of the current net node in order to detect the nearest (a) feature (designated as A), 
(b) LCS (designated as A, B, C). The aperture is a square scan containing the captured 
(a) feature A, (b) LCS ABC in its center. The aperture size is set so as to enclose at 
least one LCS. The aperture is intended for approximate determination of relative 
coordinates of the features (a) B and C, (b) A, B and C. The arrows <= between the 
features A and B, A and C symbolize skipping. With skipping, a segment is used - a 
square scan of the least sizes enclosing just one feature. 


2.1.2. Aperture scanning and LCS recognition 

At the next step, an aperture [17] scanning is carried out (see 
Fig. 2, pos. 3). The aperture is a square area around the captured fea- 
ture A which encloses several neighboring features (see Fig. 3(a)). 
During the recognition of the obtained aperture, coordinates of the 
neighboring features are determined (approximately) relative to 
feature A position. Among the detected neighboring features, fea- 
tures B and C are defined to compose, along with the feature A, the 
local calibration structure ABC |11,12,17]. 
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2.1.3. Skipping of LCS features 

After the LCS ABC has been detected, skippings (see Fig. 2, pos. 
4)|17| between feature A and its neighbors B and C are carried out, 
one after another (see Fig. 3). The skipping (designated as <) is a 
basic FOS measurement operation intended, in particular, for accu- 
rate determination of the relative coordinates of the neighboring 
features. 

The feature skipping operation A< B (skipping cycle) consists 
in moving the probe from feature A position to feature B posi- 
tion, scanning-recognizing the segment of feature B, calculating the 
“forward” differences of coordinates of features A and B, return to 
position of feature A, scanning-recognizing its segment and then 
calculating the “backward” differences of coordinates of features A 
and B. The relative coordinates of features A and B are calculated as 
a half-sum of the obtained forward and backward differences. Such 
approach permits to exclude distortions produced by the drift of the 
microscope probe relative to the sample surface [17]. 

The skipping provides accurate results under the condition of 
constant drift velocity | 17|. The skipping operation is repeated sev- 
eral times [8], and then the obtained relative coordinates of the 
features are averaged out. Once the accurate coordinates of the fea- 
tures forming LCS have been defined, the sought LCCs Ky, Ky, Kz are 
calculated (see Fig. 2, pos. 5). 

Provided that the distance between the net nodes is compa- 
rable with the distance at which the scanner is being displaced 
during the skipping, the pause insertion and the attachment-to- 
feature operation may be omitted in a number of cases. The fact is 
that the aperture scanning followed by aperture recognition is actu- 
ally a sort of attachment. Since the aperture sizes are greater than 
the sizes of the scanning area typical of a regular attachment, the 
attachment made with the help of aperture just provides a some- 
what less accurate holding of the captured feature in the center of 
the instrument’s field of view. 

For the absolute coordinates x, y, z of the LCS to which the LCCs 
Kx(X, Y, Z), Ky(X, y, Z), Kz(x, y, Z) are related, it is suitable to use the 
“gravity center” coordinates of the calibration structure. If several 
calibration structures are located near the current net node then, 
after aperture recognition, such calibration structure is selected 
among them whose coordinates are the closest of all to the node 
(see Fig. 3(a)). 


2.2. Using LCCs for determination of the true scanner movements 


Generally, the true (corrected) coordinates x, y, z of a point x, 
y, zin the image of an unknown surface are found by summing up 
the LCCs of the points of the trajectory by which the scanner has 
moved into the given image point 


Nie V5Z) = S K(x, eer 


1 


VilXj, Vis Zi) = S Ky (Xi Vis Zi), (1) 


1 


ZZ = S_KA(%. VinZ), 


1 


where Xj, y;, Z; are integer-valued coordinates of ith point of the 
movement trajectory. 

Actually, while summing up LCCs, the local unit steps corre- 
sponding to each trajectory point are summed. LCCs have a plus 
sign while moving in the positive direction and a minus sign in 
the negative direction. As in general case the LCCs depend on more 
than one coordinate, the correction result for a point of the scanner 
field is, strictly speaking, path-dependent, i.e., it would depend on 
a spatial path the scanner had passed before getting to this point. 


2.3. Types of movement trajectory by the net nodes 


In the suggested method of distributed calibration, LCCs are 
determined using the following two types of trajectory of move- 
ment by the net nodes: 


(1) With “fixed” positions of Z manipulator, coordinates of X and Y 
manipulators are changed as shown in Fig. 1(a); 

(2) With fixed positions of Y manipulator, coordinates of Z and X 
manipulators are changed as depicted in Fig. 1(b). 


In Fig. 2 (see pos. 7-9), designations at the left of a slash corre- 
spond to the trajectory of the first type, designations at the right - 
the second type. In both cases, the movements by the net nodes in 
neighboring lines/columns and in neighboring horizontal/vertical 
planes are carried out in opposite directions. Moving in opposite 
directions permits to decrease the resulting creep produced by the 
movement by the net nodes along directions x, y and Z, x, respec- 
tively. 

It should be noted that the productivity of calibration with fixed 
Y manipulator is higher in comparison with fixed Z manipulator 
since n times less aperture scans-recognitions are required for each 
net node, where n is the number of fixed positions of Z manip- 
ulator. However, as it was pointed out above, beginning with a 
certain moment, a violation of the accepted trajectory of move- 
ment by the net nodes occurs because of the long holding of the 
same feature/LCS. 


2.4. Accounting for nonorthogonality and spurious couplings 


By using lateral LCCs Kx, Ky of the orthogonal scanner [1], 
only static nonlinear raster distortions resulted from a nonlin- 
ear response of scanner piezoceramics to the applied voltage can 
be corrected. To also correct residual nonorthogonality and static 
nonlinear raster distortions resulted from local imperfections of 
movement straightness of the scanner X, Y manipulators (guidance 
errors; caused by spurious couplings Xs Y, Z> X, Z=> Y), the cal- 
ibration coefficients Kx, Ky of the nonorthogonal scanner and the 
obliquity angle w (scanner nonorthogonality) | 1-4,23-25] calcu- 
lated by LCSs in many points of the scanner field should be used. 

In the simplest case (see Fig. 3), when three features A, B, C with 
a priori precisely known distances between them are used as LCS 
[1], it is sufficient to carry out two skippings: A<~B and A<C. For 
example, it applies to the case of three neighboring carbon atoms 
on surface of HOPG monocrystal | 1 |. 

If feature sizes of the standard surface are a priori known pre- 
cisely |24], it is sufficient, instead of scanning the aperture and 
performing two skippings, to carry out counter-scanning |20| of 
the segment ofa single feature A. When making a counter-scanning, 
first, a conventional image is obtained called the direct image, after 
that a counter image is obtained by reversing the direction of the 
probe movement along a line and the direction of the probe move- 
ment from line to line |2,16,20]. The found pair of images is called 
counter-scanned images (CSIs) |20]. It was shown in Ref. [20] that 
in order to make drift correction of an SPM image possible, it is suf- 
ficient to have at least one common surface feature between the 
direct and the counter images. So, after drift in the segment has 
been corrected | 20], the true sizes of feature A may be determined, 
and then the sought for LCCs corresponding to the absolute position 
of feature A may be calculated. 

Implementation of such method of calibration in the lateral 
plane requires three different lateral sizes of the feature to be 
known |1,4]; calibration in the vertical plane requires one verti- 
cal size (for example, step height when using atoms/interstices 
as features) [2,24,26-28]. The advantage of calibration by known 
feature sizes consists in a substantially greater productivity. 
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Provided that both the distance between LCS features and the sizes 
of the features of the standard surface are all known, the precision 
of distributed calibration may be improved. To do this, instead of a 
regular scanning during the skipping operation, a counter-scanning 
of the segments should be performed followed by determining the 
LCCs by the above two methods. 

The distribution of LCCs in the plain field is obtained by perform- 
ing scanner calibration for some “fixed” position of Z manipulator. 
Beside nonlinearity and nonorthogonality of the scanner, mutual 
spurious couplings of X<Y type are also taken into account in 
the found distribution. The change in Z manipulator position while 
moving in the lateral plane from one LCS to another shows the 
value of error induced by spurious couplings of X > Z, Y=>Z types 
(the bowl-shaped image effect well noticeable on flat samples) 
[11,24,25,29]. 

The spatial distribution of LCCs Kx, Ky can be obtained by per- 
forming scanner calibration for different positions of Z manipulator 
set by the coarse approach stage. Beside nonlinearity in the ver- 
tical plane, the spurious couplings of Z>X and Z=>Y types will 
be accounted in the found distribution. In general, the spurious 
couplings Xs Y, Zs=X, Z<Y arise due to imperfections in con- 
struction, control, materials, and manufacturing process of the 
scanner. 


2.5. Calibration database 


The LCCs obtained during the distributed calibration along with 
the absolute coordinates of LCS to which they correspond are stored 
in the calibration database (CDB, see Fig. 2, pos. 6). By storing the 
dates of the local calibrations in a database, it is possible to accu- 
mulate the information on repeated calibrations for several years, 
which allows determining the scanner piezoceramics age dynam- 
ics |8|. As the age dynamics are known, it is easy to find such time 
interval during which the accumulated LCCs would yield a regres- 
sion surface (see Section 3) ensuring correction of the surface image 
with some minimum error. 


2.6. Adaptive properties of the algorithm 


The modulus of drift velocity during attachment, skipping or 
counter-scanning should not exceed a certain value — the one at 
which the current feature can go out of the limits of the scanned 
segment. Another condition imposed on the drift is invariability of 
the drift velocity during one cycle of skipping or segment counter- 
scanning | 17,20]. The better the last condition is satisfied, the less 
error of drift elimination may be obtained which, in its turn, would 
lead to a more accurate LCC determination. Measurements show 
that the drift velocity remains practically constant over tens of 
seconds | 16,17,20,21] or even tens of minutes | 30-32] whereas one 
cycle of skipping between neighboring carbon atoms on pyrolytic 
graphite surface may take as little as 300 ms | 17]. 

To eliminate the calibration errors connected with large changes 
in drift velocity, a ceaseless in situ monitoring of drift velocity 
is established while skipping and/or counter-scanning |17-21]. 
Once an unallowable change in drift velocity is registered, the data 
obtained during the current skipping cycle are declared spoiled and 
discarded. After that, several idle skipping cycles are inserted. If 
the drift velocity has recovered during that time, the interrupted 
skipping is continued, otherwise it is restarted. In case of a long- 
term instability, the microscope probe is returned to the current 
net position and the local calibration is carried out once again. 

If an unallowable change of the drift velocity interrupts the 
calibration by feature sizes, then an idle counter-scanning of the 
segment is applied. Here, like with the skipping, the local calibra- 
tion for the current net node will be restarted if the drift velocity 
is unable to stabilize during the set time interval. The calibration 


algorithm responds the same way to an increase in noise, which 
may derive from one of the following: surface defect, contami- 
nation of tip/surface, activation of a vibration source, change of 
humidity, electromagnetic disturbance, etc. The suggested calibra- 
tion algorithm can thus be considered as adaptive, as it changes its 
operation in accordance with the actual changes in measurement 
conditions. 


3. Correction of scan of an arbitrary surface 


Because of the drift as well as the differences in dimension, 
structure, and orientation between the cells of the initial net and the 
lattice of the standard surface, the obtained LCC net is no longer a 
square cell net with integer period, and so it should be transformed 
to one for subsequent use. The transformation consists in construc- 
ting a regression surface, i.e., drawing through the noisy LCCs such 
a smooth surface f that a sum of squared deviations between the 
LCC values and the corresponding surface points is minimal. These 
constructions may be carried out, for instance, for each fixed posi- 
tion of the scanner Z manipulator. As a result, the following set of 
surfaces is obtained 


Ki = fr, (x, y, Zi), 
Ry = fi (%Y. 2 
a’ = fa(x, y, Zi), 

Kz = fx,(%, Y, Zi), 


where letter r means regression and Zz; are fixed positions of the 
scanner Z manipulator (i=1,2,...,n) adjusted during calibration by 
means of coarse approach stage displacement. Constructing regres- 
sion surfaces (2) also permits to reduce the influence of errors of 
determination of LCCs and local obliquity angles | 1] on the results 
of the nonlinear correction. 

Formally, in any (x, y) point the pair of coefficients Ky, Ky of some 
hypothetical orthogonal scanner may be expressed in terms of the 
calibration coefficients Kx, Ky, and the obliquity angle aw by using the 
following transformation that establishes a relationship between 
the unit length scales of a rectangular and an oblique system of 
coordinates | 1 | 


Kx(x, y) = K(x, y) + Ky(x, y) sin [a(x, y)], a 
Ky(x, y) = Ky(x, y) cos [a(x, y)]. 


The z coordinate in transformations (3) is omitted to simplify the 
notation. 

By replacing LCCs Kx, Ky, and obliquity angle w in formulae (3) 
with the values taken from the corresponding regression surfaces 
(2) and then substituting the obtained expressions for Kx, Ky in (1), 
one may write 


Xi(Xj, Vi) = SX ¥i) + Ky (xj, yi) sin [a" (xj, Yi) I, 


a (4) 
YilXi, Yi) = S Kj (%i, yi) cos [a"(x;, vi), 


1 


where X;, y; are integer-valued lateral coordinates of the ith point 
of the movement trajectory. 

Strictly speaking, the K"(x;, y;) term in transformations (4) is 
responsible for a shift of the point to be corrected along x if only 
the movement trajectory of the scanner has an x-component at this 
point. The Kt (x;, y;) sin [a@'(x;, y;)| term sets a shift of the point to be 
corrected along x if only the movement trajectory of the scanner 
has an y-component at this point. The Ky (X;, ¥;) Cos [a"(x;, y;)| term 
provides a shift of the point under correction along y if only the 
movement trajectory of the scanner has an y-component at this 
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point. Actually, the above conditions establish the rules for cal- 
culating the corrected coordinates x, y. To represent these rules 
analytically, formulae (4) should be rewritten as follows 


(xi Vi) = > R(x YAK + Kj (x, yi) sin la” (x, yi)]Ayi. 


_ (5) 
Hilxi Vi) = Y_ K(X, Yi) cos [0 (xi, YAY. 

1 
where Ax; =X; —X;_1, AY;=¥;—Yi_1 are differences of integer- 
valued coordinates of the trajectory neighboring points, which may 
take the values O or +1. 

It might seem that we could use one regression surface less if we 
calculate LCCs Kx(x, y), Ky(x, y) according to formulae (3), replace 
them with corresponding regression surfaces Ky(x, y), K(x, y), and 
then substitute in (1). However, we must not do this. The fact is 
that the “switches” Ax;, Ay; provide the information of the trajec- 
tory, which would be fixed in this case. So the CDB will be bound 
to a single particular trajectory and, as a result, will lose its uni- 
versality and therefore its practical value. Thus, LCCs Ky, Ky, and 
local obliquity angle a should be stored in the CDB and processed 
separately. 

In formulae (5), polynomials are suitable for use as regression 
surfaces 


Ki (x, Y) = do + G1y + d2X + G3xXy + d4y? + a5x? +--- =~ axPax yPay, 
K!(x, y) = bo + Diy + box + b3xy + bay? + b5x2+--- =S~ bxPbxyPoy, 
al" (X,Y) = Co + Cy + CoX + C3XY + Cay? + C5X27 +--+» =)_ exPexyPey, 
KI (x, y) = do + diy + dox + d3xy + day? + d5x*+--- => ar 
6 


where a=(do, 41, d2,...),D=(bo, by, b2,.. .), C=(Co, C1, C2,...), d=(do, 
d,,d,...) are vectors of the polynomial coefficients; px =(0, 0, 1, 1, 

2,...), Py =(0, 1, 0, 1, 2, 0, ...) are vectors of powers of x, y vari- 
ables, respectively. The highest degree of the polynomial for each 


Formulae (7) are intended for correction of raster trajectory 
of the direct image (designated as number 1 in Fig. 1(a) in [20]). 
Proceeding by analogy, it is easy to obtain formulae that correct 
the raster trajectory of the counter image (designated as number 
2 in Fig. 1(a) in |20]). To do that, one should consider the opposite 
direction of the trajectory and its origin location in the coincidence 
point (Xmax, Ymax) [20]. 

With raster scanning, the trajectory pattern is known a priori 
and may be easily formalized. As retrace is used in a line, it allows 
passing from trajectory points x;, y; to arbitrary raster points x, y, 
the formulae (7) take the form 


X(x,y) = SOKi(x0 +X, Yo +) 
x 
+S Kyo, Yo +y)sin[a" (xo, Yo + YI. (3) 
y 

Wy) =) _Ky(X0, Yo +) cos [a" (Xo, Yo + Y)]- 

y 
Here, Ax, Ay differences are present in formulae (8) implicitly 
through the known movement trajectory in the raster. 

Instead of sum accumulation in expressions (8), the sought coor- 


dinates may be determined directly by calculating the following 
integrals 


x 
X(x,y) = / Ky(Xo +X, Yo + y)dx 
0 
y —_ 
a / Ky (Xo, Yo + y) sin[@" (Xo, Yo + y)]dy, (9) 
0 
y 
yy) = / Ky (Xo, Yo + y) cos [a@"(X0, Yo + y) | dy. 
0 
The formulae for counter image are obtained likewise. 
For regression surfaces represented by polynomials (6), the inte- 


grals (9) have analytical solutions. For planes, for example, the 
following formulae can be obtained 


: 1 
X(X, Y) = [do + 41(Yo + Y) + G2Xo |X + ox + — = |2 sin ( a) { (bo + bio + b2X9)sin ec + Cy (vo + ) + €2X0 | 


by 


1 


regression surface (6) may be chosen so as to ensure, for exam- 
ple, the least error of the mean value of the calibration size on the 
corrected image of a standard | 14]. 

Transformations (5) are universal and applicable for coordinate 
correction of any kind of trajectory. When the trajectory is a raster, 
transformations (5) take the following form 


Xi(Xi, Vi) = S (K(X + Xi, Yo + Vi JAX; 
+K}(Xo, Yo + yi) sin[a" (Xo, Yo + Vi) Ayi}, (7) 
YilVi) = SK} (Xo, Yo + y;) cos [a"(X9, Yo + Vi)IAYi, 


1 
where Xo, Yo are Coordinates setting araster position in CDB relative 
to the scanner origin of coordinates; x;, y; are points of the raster 
trajectory of the distorted image which belong to the ranges 0, 1, 2, 
..Xmax and 0, 1, 2,..., Vmax, respectively (Xmax + 1 is the number of 
points in the raster line and ymax + 1 is the number of raster lines). 


+ <, 60S cc +Cy (vo + a + 2X0 — byycos[co + €1(¥Vo +y) + co] : 
1 


2 1 
y(y) = Gi 2 sin ( oY) { (bo oe bio ae b2Xo9 )cos cc + Cy (vo ae ) f C20 


(10) 


bi. | 
~ sin ec +1 (v0 + a + 2X0 \ + biysin[co + c1(yo +y) + 2X0 


For higher order polynomials, the analytical solutions become 
yet more cumbersome. To avoid intricate formulae, the inte- 
grals can be calculated by applying numerical methods. Moreover, 
by using in the integrals (9) instead of Ki (x, y)sin[a’™(x, y)] and 
oe y) cos [a'(x, y)| the corresponding regression surfaces 


(X,Y) = > exPex yPey | 
(X,Y) = >> fxPex yPty 


joe e, f and p are defined similar to those in (6)) built respec- 
tively for products Ky(x, y) sin [a(x, y)] and Ky(x, y) cos [a(x, y)] (the 
analogous substitution may be used in formulae (7), (8)), a more 
compact formulae may be written down 


(11) 


. a 
en ne x)Pat! — xoPax*l](yo + y)Pav 


as 





xoPex[(yo + y)Pey*! — yoPey*"], (12) 


Pey + 1 





XoP&[(Yo + y)POt! — yoPtyt!]. 


Pry + 1 
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Strictly speaking, the numerical values obtained by formulae (12) 
are not equivalent to the values obtained by formulae (9) but are 
pretty close to them. 


4. Estimation of scanner errors 


The more symmetrical scanner construction is used, and the 
more precise its fabrication, more homogeneous material, and 
more evenly distributed load, the more symmetrical the surfaces 
(2) are and the stronger they resemble solids of revolution. Asym- 
metry or deviation from a revolution solid may serve as signs of 
constructional, technological or material imperfection of the scan- 
ner or incorrect scanner control. In absence of nonlinearities and 
spurious crosstalk couplings, the regression surfaces (2) degenerate 
into horizontal planes. Thus, the difference between the obtained 
regression surface and the horizontal plane may serve as a measure 
of the nonlinear distortions and spurious couplings of the scanner 
[8,14]. To get estimate for this measure, the following maximum 
differences will be calculated 


Amax — max |K7(x, y) — K}(0, 0) 





’ 


Amax — max K(x, y) —K7(0, 0)|, (13) 


kK" 
y 
Amax — max Jar" (x, y) — a"(0, 0) 








? 


where point (0,0) corresponds to the origin of coordinates of the 
direct/counter image in case of the virtual calibration | 14] or to the 
scanner’s origin of coordinates in case of the real calibration | 8]. 

To estimate root-mean-square deviations of LCCs and obliquity 
angle the following formulae will be used 








n 

1 > - 2 

OK, = > [Kelis Yi — Kei Yi] ; 
i=1 
1 n 

- : 2 

OG = [Ky (xi. vai — KX Yi) (14) 

i=1 





n 
1 
ar | ) [a(X;, Vi)i — a(x, vil, 
i=1 


where Xj, Vi, Kx(X;, Vii, Ky(%, Vi); @(Xj, Yi); are the coordinates, LCCs, 
and the local obliquity angle of the ith LCS, respectively, taken from 
the CDB; n is the number of LCSs in CDB. 


5. Application of feature-oriented positioning technique 


In a number of practical cases it is quite enough to calibrate 
precisely just a small part of the scanner space located in the vicinity 
of origin of coordinates, where the residual errors are minimal [8]. 
To perform the required measurements, the area of interest on the 
surface can be simply transferred to this part of the scanner space 
by applying feature-oriented positioning (FOP) methods | 17-19]. 

A successful solution of the problem of finding LCCs distributed 
in movement space of the scanner mostly depends on the prop- 
erties of the standard surface used. The latter should possess an 
invariable structure in each point of the scanner field. In practice, 
however, defects and residual stresses distort the surface structure 
of the standard somewhat complicating the distributed calibration 
on large fields. In this case, a preliminarily selected “perfect” region 
of the standard surface has to be moved within the scanner space by 
using the same FOP methods in order to calibrate the whole space 
only by this small area of the standard surface | 17-19]. 

In both cases to perform the transport operations mentioned 
above, a coarse X, Y, Z positioner is required capable of carrying 


out a long-distance travel with the step comparable to the typical 
lateral size of the used reference feature | 17-19]. It is important 
to note that no special requirements of movement precision are to 
be met by the coarse positioner, since with FOP the entire move- 
ment precision is provided by the scanner (fine positioner) and its 
calibration | 17-19]. 


6. Discussion 


Yet another variation of distributed calibration is shown in 
Fig. 3(b). Here, during the probe attachment, the whole LCS is 
involved rather than a separate feature, and the aperture contains 
one LCS only. In fact, among the features recognized in a local scan, 
those are being searched for which make up the LCS nearest to the 
current net node; after that an attachment of the local scan is car- 
ried out to the gravity center of the found LCS. This variation of 
the distributed calibration has a higher productivity since it uses a 
reduced-size aperture. The method is applicable when the step of 
the initial net is greater than a typical size of LCS. 

In case the step of the net is comparable with the typical LCS 
size (see Fig. 3(a)), it is possible after skippings A B and A<=C to 
carry out additional skippings between feature A and all the other 
features detected in the aperture and forming an LCS with feature 
A. The obtained LCCs should be associated with the coordinates of 
the gravity centers of the corresponding LCSs. For instance, accord- 
ing to the configuration shown in Fig. 3(a), 6 times less number 
of movements between the net nodes, waitings, attachments and 
aperture scannings may be executed during that kind of calibration. 

The skipping operation may be omitted at all when no high pre- 
cision calibration of the scanner is required yet a sizeable scanner 
field is to be calibrated or when a less accurate calibration is to 
be performed for a shorter period of time. By applying counter- 
scanning in the aperture |20| instead of skipping, it is possible to 
correct drift in the aperture and then determine the sought for LCCs. 

With the calibration method suggested (see Fig. 2), the number 
of skipping cycles (pos. 4, see also Ref. [8]) are not recommended to 
be set large. Otherwise, the microscope probe may shift at a large 
distance off the current net node due to drift that after moving the 
probe to the next net node will induce a strong creep resulting in 
a long pause. Here, the trajectory of movement by the net nodes 
accepted in this method (see Fig. 1) is violated which generally 
leads to undesirable “runaway” of the scanner’s piezomanipulators. 
Moreover, with the large number of skipping cycles, the relative 
distances between the features measured in a sequence of repeated 
cycles will correspond, owing to drift, to different absolute pos- 
itions of the scanner, where, strictly speaking, different nonlinear 
Static distortions occur. 

Thus, to obtain more accurate values of the noise-sensitive 
calibration coefficients Ky, Ky and the obliquity angle a [1], the cali- 
bration routine should be repeated many times (see Fig. 2, pos. 10) 
adding each time a new set of LCCs to the database. The same is true 
with respect to a number of consequently executed counter scans 
obtained calibrating by sizes of feature of the standard surface. 

Similar arguments exist against the attempts to measure LCS 
with redundancy by performing an additional skipping B < C fol- 
lowed by least-squares data adjustment as it is accepted in geodesy 
[33], for example. Nevertheless, for a small drift, the use of this 
trick allows to increase the precision of measurement of feature 
coordinates. 


7. Conclusion 


The developed calibration method employs a number of prin- 
ciples and tricks embeded in the FOS methodology |17], namely: 
handling separate surface features, movement by short distances 
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from one feature to another located nearby, measurement of the 
relative distances, ultimate localization of all measurements, mul- 
tiple repetitions of the measurements, ceaseless probe attachments 
to the surface features, continuous monitoring of the drift velocity, 
drift distortion neutralization by means of hierarchically-organized 
counter movements. Thus, the suggested calibration method can be 
thought of as another component of the FOS methodology. 
The key points of the research can be summarized as follows: 


(1) A method of automatic distributed calibration has been 
suggested that allows excluding errors caused by static nonlin- 
earities, nonorthogonalities and spurious couplings of a probe 
microscope piezoscanner. 

(2) The calibration database obtained during calibration is insensi- 
tive to microscope head thermal drift and piezoscanner creep 
due to the feature-oriented scanning methodology applied. 

(3) The lattice constant of a crystal has been used as a nature length 
standard of nanometer range. 
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